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ABSTRACT 



Aims. We study the properties of X-ray emitting plasma of MP Mus, an old classical T Tauri star. We aim at checking whether an 
accretion process produces the observed X-ray emission and at deriving the accretion parameters and the characteristics of the shock- 
heated plasma. We compare the properties of MP Mus with those of younger classical T Tauri stars to test whether age is related to 
the properties of the X-ray emission plasma. 

Methods. XMM-Newton X-ray spectra allows us to measure plasma temperatures, abundances, and electron density. In particular the 
density of cool plasma probes whether X-ray emission is produced by plasma heated in the accretion process. 

Results. X-ray emission from MP Mus originates from high density cool plasma but a hot flaring component is also present, suggesting 
that both coronal magnetic activity and accretion contribute to the observed X-ray emission. We find a Ne/0 ratio similar to that 
observed in the much younger classical T Tauri star BP Tau. From the soft part of the X-ray emission, mostly produced by plasma 
heated in the accretion shock, we derive a mass accretion rate of 5 x 10"" Mq yr"'. 
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1. Introduction 
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Low mass stars are sources of strong X-ray radiation since 
their early evolutionary phases. Coronal plasma, responsi- 
ble for the X-ray emission, is confined and probably heated 
by magnetic fields which emerge from the stel lar surface 
dFeigelson & Montmerle"1999'; 'Preibi sch et al.ll2005l) . The coro- 
nal plasma observed in essentially all late-type stars can be char- 
acterized by a large variety of average temperatures (from few 
MK to tens of MK) and metallicities (from one tenth to few times 
the solar photospheric value), but a common feature is the low 
density measured at the temperature of formation o f O vii He- 
like triplets (A^e _~ 10' " - 10" cm"^ atT ~ 2MK) jTesta et all 
l2004tlNess et"ani2004 . 

In very young stars, however, accretion may cause X-ray 
emission in addition to magnetically-confined coronal plasma. 
In classical T Tauri stars (CTTSs) gas falls from the circumstel- 
lar envelope, funnelled by the magnetic field, and hits the stel- 
lar photosphere. In the resulting sho ck the accreted material i s 
heated to temperatures of few MK ( Calvet & G ullbring' ll998l) . 
With typical mass accretion rates, the shock-heated plasma can 
reach X-ray luminosities as high as 10-" erg s"'. 

The few CTTSs for which high-resolution X-ray spec- 
troscopy was performed up to date show, in most cases, cool 
plasma components (2 - 4MK) with large electron densi- 
ties (10" - 10'^ cm""*) which have been interpreted as evi- 
dence for X-ray emission due to an accretion shock. The best 
known examples of th is behavior are the CTTSs TW Hya, 
BP Tau, and V4046 Sgr (iKastner et al.l2002HSchmitt et al.l2005t 
iGiinther et al.ll2006l) . Noticeable exceptions are the CTTS T Tau 



riudel et alJl2007h and the Herbig star AB Aur dTelleschi et al.l 
I2007h . 

In the hypothesis of X-ray emission originated in shocks 
great attention has been focused also on plasma element abun- 
dances. In fact, they probe the chemical composition of the ac- 
creting stream, and hence provide insightful indications on the 
physical and chemical processes at work in the inne r circumstel- 
lar disk ( Stelzer & Schmittll2004l:rDrake et aljllool) . 

In this letter we present the XMM-Newton observation of 
MP Muscae, one of the oldest known CTTSs, aimed at studying 
the properties of the X-ray emitting plasma and the role of the ac- 
cretion process. MP Mu s is a Kl IVe star of the Lower Centaurus 
Crux (LCC) association. iGregorio-Hetem et al. I (ll992h first iden- 
tified it as a classical T Tauri star by measuring enhanced Ha 
emission (EW = -47 A) and Li absorption (EW = 0.37 A). 
Excesses in infrared bands revealed an optically thick cir cum- 
stellar disk dMamaieket all 120021; ISilverstone et al.ll2006l) wit h 
an estimated dust mass of ~ 5 X 10-^ M© ( Carpenter et al.l2005h . 
Batalha et al. ( 1998) derived a rotational period of 5.75 d from 
variability in the B, V, R, and I bands, although the amplitude of 
var iability is surpri si ngly l ow for a CTTS. 

iMamaiek et alj (l2002h derived three different ages for 
MP Mus, 7, 14, and 17Myr, depending on the adopted the- 
oretical evolutionary tracks. The LCC association, to which 
MP Mus belongs, is the oldest portion of the Scorpius-Centaurus 
OB association, and its estimated age is b etween 16 and 23Myr 
dMamajek et aI.ll2002l:[Sa-tori et al.ir2003h . 



2. Observation and data analysis 



Send offprint requests to: C. Argiroffi, e-mail: argi@astropa.unipa.it 



MP Mus was observed with XMM-Newton for a duration of 
~ llOks on 2006 August 19-20. We processed the data using the 
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Fig. 1. Background-subtracted light curve of MP Mus obtained 
by adding the three EPIC instruments. 

Table 1. MP Mus best-fit paramenters. 
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Fig. 2. RGSl spectrum in the wavelength region of the Ovii 
triplet (gray) with the best fit lorentzian line profile (black). 
Wavelength bins corresponding to bad column pixels are plotted 
without error bars and were not considered in the fitting proce- 
dure. 
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to the 68% confidence level. 



SAS V6.5 standard tasks. After having discarded time segments 
affected by high background count rates, we obtained good time 
intervals summing up to ~ 100 ks. The X-ray light curve (Fig.[TJ 
shows clear hints of flaring activity, typical of magnetically ac- 
tive coronal sources. 

To increase the signal-to-noise ratio we rebinned the PN and 
MOS spectra to obtain at least 30 counts in each bin, and the 
RGS wavelength bins were joined two by two. Spectral anal- 
ysis was performed usin g the Astrophysica l Plasma Emission 
Database (APED V 1.3, fSmith et alj lioOll) . We adopted the 
lAsplund et al.l ( l2005h solar photospheric composition as the ref- 
erence set of element abundances. 

We derived the characteristics of the X-ray emitting plasma 
of MP Mus by fittin g PN and MOS spectra with XSPEC 
VI 1.3.0 (lArnaudll996h . adopting an optically-thin plasma emis- 
sion model with three isothermal components. The best-fit 3-T 
model (Table [1]) also provided individual O, Ne, Fe, and S abun- 
dances, while the abundances of the other elements were tied to 
the Fe abundance because the fit quality does not improve if they 
are treated as additional free parameters. 

Table |2] shows the fluxes of the strongest R GS lines, mea- 
sured with PINTofALE V2.0 dKashvap & Dra ke 2000), adopt- 
ing a lorentzian function to fit the line profile. We checked that 
the derived 3-7" EPIC model provides a reasonably good descrip- 
tion of the RGS line fluxes. 

We identified the O vii and Ne ix He-like triplets in the RGS 
spectra. The O vii lines (indicated as r, /, and / in Fig.|2]i provide 
a density-sensitive ratio /// = 0.28 + 0.13, which implies an 
electron density log A^e - 11-7 + 0.2 for the plasma at T ~ 2 MK. 
The Ne ix line flux measurements are affected by large uncer- 
tainties due to the strong blending with Fe lines. Figure [3] shows 
the observed Ne ix triplet with that predicted on the basis of the 
3-T EPIC model assuming different A^e values. This comparison 
suggests that the plasma at T ~ 4 MK, which produces the Ne ix 
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Fig. 3. RGS2 spectrum around the Ne ix triplet (gray) with the 
predicted spectra obtained from the 3-7" EPIC model and assum- 
ing different electron densities. 



triplet, has a density logA^e 
estimated from the O vii. 



3. Discussion 



> 11, and possibly as high as that 



From the analysis of the O vii triplet we find that the cool plasma 
component in our source has a density significantly larger than 
typical coronal values. It suggests that shock-heated plasma con- 
tributes significantly to the observed X-ray emission. In this re- 
spect, MP Mus is the fourth CTTS discovered so far with evi- 
dence of X-ray emission produced by cool high density plasma, 
likely resulting from an accretion process. Previous cases were 
TW Hya, BP Tau, and V4046 Sgr Moreover, MP Mus shows 
also clear evidence of intense coronal activity, as indicated by 
the flares (see Fig.lTJ and by the hot plasma component. 

3.1. Abundances 

The X-ray emitting plasma in MP Mus is heavily depleted of 
Fe, Si, and Mg; O, and S are moderately depleted, while Ne dis- 
plays a larger abundance. If we assume that soft X-rays are pro- 
duced by plasma heated in the accretion process, the observed 
abundances probe the chemical composition of the infalling cir- 
cumstellar material. We compared the present abundance val- 
ues of MP Mus with those obtained for the other three CTTSs 
showing evidence of X-ray emission due to shock-heated plasma 
(Table|3]l. In all cases the X-ray spectra indicate that the accreted 
material has a Ne abundance enhanced with respect to the other 
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Table 2. Strongest RGS lines of MP Mus. 



Table 3. CTTSs properties. 
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Line fluxes (10 phs cm -) with uncertainties at 



elements (with respect to the solar photospheric abundance ra- 
tios); 

IStelzer & SchmittI ( |2004|) suggested that the non-solar abun- 
dances of the shock-heated plasma of TW Hya may be explained 
by assuming that the accreting materia l underwent grain dep le- 
tion, a mechanism already invoked by iHerczeg et alj (l2002h to 
explain the low Si abundance. Depending on the temperature, the 
circumstellar material is composed by gas and dust phases which 
have different chemical compositions. The actual abundance ra- 
tios in the gas phase are determined by the di fferent condensation 
tempe ratures of the various elements dSavage & SembachI 
Il996l) . One of the proposed scenarios of disk stratification pre- 
dicts that dust grains mostly settle in the disk midplane, while 
the gas extends up to the disk surface, where it is easily ion- 
ized by the stellar radiation. Here it is funnelled along the mag- 
netic field lines and accretes onto the central star. If the gas and 
dust phases are spatially separated in the inner circumstellar disk 
and/or the gas accretes more efficiently than the dust, the shock- 
heated material should be depleted of those elements (like Fe) 
which easily condense into dust grains and conversely enriched 
of more volatile elements (like noble gases). Following this sce- 
nario, the accreted material should have an abundance pattern 
similar to that observed in the interstellar gas, i.e. with the abun- 
dances decreasing for increasing T^. Note that the phenomenon 
of gas/dust separation in circumstellar disks, and the subsequent 
accretion of only the gas phase, is also hel d responsible for the 
peculiar abundances of RV Tauri stars (see lGiridhar et alJl2005L 
an d references therein ) . 

iDrake et all ( l2005h showed that TW Hya displays a Ne/O ra- 
tio larger by a factor ~ 2 than the uniform value observed in 
a vast sample of coronal X-ray sources (Drake & Testa 2005). 
Moreover, Ne/Fe abundance ratios as large as that of TW Hya 
were observed only in very few cases of putatively coronal 
sources. These results support the hypothesis that soft X-ray 
radiation from TW Hya is not produced by coronal plasma. 
[Drake et al. suggest that the Ne/O ratio can be used as a cri- 
terion to identify X-ray sources where the emitting material in 
the accretion stream suffered grain depletion. However, we note 
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Fig. 4. Average plasma temperature and Ne/O ratio vs age for 
the sample of four CTTSs with evidence of high density cool 
plasma. 



that the condensation temperature of O is quite low (180K), 
therefore the separation between gas and dust must occur at low 
temperature to produce significant O depletion in the accretion 
streams. 

A large Ne/O abunda nce ratio is observed also in the X-ray 
spectrum of V4046 Sgr (iGtinther et al.ll200 6'). where high den- 
sity hints again at X-rays from shock-heated plasma. Instead, 
both MP Mus and BP Tau have Ne/O ratios typical of stellar 
coronae. Drake et aTl ( l2005h explained the Ne/O ratio of BP Tau, 
lower than that of TW Hya, on the basis of the different evolu- 
tionary stages of their circumstellar disks. Since TW Hya is sig- 
nificantly older than BP Tau, it is conceivable that the dust/gas 
separation process, and the subsequent depletion of high el- 
ements, is not visible in the latter case because these processes 
occur on a time scale longer than the age of BP Tau (~ 0.6 Myr). 

In Table [3] we report the ages of the four CT TSs introduced 
above. We adopt an age of 17 Myr for MP Mus (iMamaiek et alJ 
I2OO2I) . since it is compatible with the age of the LCC association. 
For the subsequent discussion, the absolute age of each CTTS is 
unimportant, while only the age sequence matters, whose reli- 
ability depends only on the correctness of the membership of 
these CTTSs to the relevant stellar associations. 
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Figure |4] shows the variations of plasma average tempera- 
ture and Ne/O ratio with respect to stellar age, for the sample 
of four CTTSs (having spectral types ranging from Kl to K7). 
Both Tmed and Ne/O do not have a monotonic trend with age, 
but these two plots suggest that stars with hotter plasma have 
lower Ne/O ratios, and vice versa. It is likely that high r^ed in- 
dicates a large contribution from coronal plasma to the whole 
X-ray emission. In this scenario of mixed accretion-driven and 
coronal X-ray emission, the measured Ne/O ratio is a weighted 
average of the values in the shock-heated plasma and in the coro- 
nal plasma. Hence any large Ne/O ratio of the accreted mate- 
rial may be partly hidden by the coronal plasma abundances. 
To check this possibility we fitted the observed EPIC spectra of 
MP Mus assuming a high Ne/O ratio for the coolest plasma com- 
ponent, but the model does not reproduce the observed spectra 
as well as the model described in Sect. |2l Moreover, in MP Mus 
the hot coronal plasma does not contribute significantly to the 
observed O and Ne line emission (see below). 

We conclude that the relatively low Ne/O ratio in MP Mus 
is a characteristic of the cool accretion component, and the stel- 
lar age is likely not the only parameter which determines the 
Ne/O ratio observed in CTTSs with evidence of high density 
cool plasma. 

3.2. Accretion 

For the subsequent discussion we first assume that the cool X- 
ray emitting plasma of MP Mus is only due to the shock accre- 
tio n, with no contr ibution from coronal plasma. Starting from 
thelM amaiek et al.' (2002) results on MP Mus, and based on the 
ISiess et al.. (.2000) stellar models, we adopt for MP Mus a mass 
of 1 .2 M0 and a radius of 1 .3 Rg- 

Using the O vii triplet and the O viii Lya lines we infer the 
electron density (A^e = 5 x 10" cm""*), temperature (T - 3MK, 
obtained from the Oviii Lya and Ovii r lines), and emission 
measure (EM = 2.4 x 10^"^ cm"^) of the post shock plasma. In 
the strong shock scenario, the relevant plasma parameters are 
linked by the relations: 

1 3 umu 

Ni = 4No, VI = -vo, Ti = 77^v^ (1) 
4 LO k 

where the suflixes and 1 indicate the pre-shock and post-shock 
plasma, the density, v the velocity, T the temperature, and fi 
the mean molecular weight (in our case fj. - 0.61). From the 
measured temperature T[ we infer that the pre-shock velocity 
is 470 km s"'. This value corresponds to a free fall from an in- 
ner radius of the circumstellar disk of 3/?*, or from a larger 
distance if some energy loss occurs during the fall. From the 
post-shock plasma temperature and density we derive a cool- 
ing time of 350 s, and considering that the post-shock velocity is 
120 km s"', we obtain a characteristic length of the post-shock 
region Z = 4 x 10^ cm - 0.05 Ri,. Hence, the cross section of the 
infalling stream A - EMI{N^N\il) is 3 x 10^" cm^. It corresponds 
to a filling factor / = A/(4;r/?^) of 0.3 % of the stellar surface, 
and to a mass accretion rate of 5 x 10"" Mq yr"^ 

We made the hypotheses that: (1) the cool plasma is pro- 
duced in the accretion shock; (2) the cool plasma is optically 
thin; (3) its density is measured from the O vii ///. 

We are confident that the assumption (1) is appropriate. First 
note that the two strong flares detected, produced by coronal 
plasma, contribute just 3.6% of the spectrum above 18 A (i.e. 
below 0.7 keV); second, we tried to fit the Ovii triplet with 
two contributions due to low and high electron density (10^ and 



10'^ cm respectively), finding that at least 80% of the O vii is 
due to high density (i.e. shock-heated) plasma, and at most 20% 
to low density (i.e. coronal) plasma. 

Hence the derived accretion rate M, which depends only on 
the hypothesis (1) and (2), but not on A^e, is acceptable, but a 
larger M could be possible if part of the X-ray emission is ab- 
sorbed. 

The measured A^e is more uncertain: a small contribution of 
low A^e coronal plasma to the O vii triplet might cause an under- 
estimation of A^e; conversely an UV field might influence the 
populations of the Ovii atomic levels by photoexcitation and 
hence mimic an high density plasma. No UV excess emission, 
which could originate from the accretion hot spot, has been re- 
ported for this star. However a sufficiently high UV radiation 
density can be present only very near the accretion hot spot 
on the stellar surface, and the photoexcitation hypothesis would 
anyway indicate that the cool X-ray emitting plasma is close to 
the base of the accretion funnel. 



4. Conclusions 

From the analysis of the XMM-Newton observation of the CTTS 
MP Mus we derived evidences that plasma heated in the accre- 
tion shock produces the soft part of the X-ray emission. We mea- 
sured a Ne/O ratio similar to that of BP Tau and reduced by a 
factor 2 with respect to that of TW Hya and V4046 Sgr: this 
result suggests that the stellar age is not a useful parameter to 
predict the amount of grain depletion sufltered by the accreting 
material. 
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